Eigenvalue correlations on Hyperelliptic Riemann surfaces. 
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Abstract 

In this note we compute the functional derivative of the induced charge density, on a thin 
■ conductor, consisting of the union of g + 1 disjoint intervals, J := U|j^(aj, bj), with respect to 

an external potential. In the context of random matrix theory this object gives the eigenvalue 
fluctuations of Hermitian random matrix ensembles where the eigenvalue density is supported 
(SJ ' on J. 
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1 Introduction 



Consider the minimization problem 



(1) 



E 



inf 



log \x 



where the set A consists of all positive Lebesgue measures n(t)dt such that J fi(t)dt = 1. The above 
formula describes the electrostatic equilibrium in which charges are placed on the real line in the 
presence of an external field v(x). For analytic external fields it is well known that the infimum 
(|l|) is attained at a unique measure a(x)dx which is called the induced charge density. Moreover, 
the support of the induced charge density is generically characterised by a finite number of disjoint 



intervals J := U?t^(a,j,bj] 
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At electrostatic equilibrium, the induced charged density a(x), x € J, satisfies the following 
integral equation 



(2) v(x) — 2 J In \x — t\a(t)dt = A = constant, x € J. 

where the constant A is the Lagrange multiplier which fixes the constraint 

(3) J a(x)dx = 1. 

To determine a, we convert (2) into a singular integral equation by taking the derivative of @ 
w.r.t. x, that is 

(4) 2P [ ^ = x e J, 

J j x — t ax 

where P denotes the principal value of the singular integral. For generic potential v such that v'(x) 
is Holder continuous]] the solution, cr, of the singular integral equation (|j), which is bounded at the 
end points of J, is necessarily zero there |nj , i.e. 

a(aj) = = a(bj), j = l,...,g + l. 

The end points {oj, of the support of a are determined by @, 

(5) / a{x)dx = 0, j = 1,2,..., g, 
Jb 3 

and by the moment conditions, 

(6) / dx = 0, k = 0,...,g. 



See P]. The equation (|2|) also arises from a mean-field approach to random Hermitian matrix 
ensembles §; a(x), is the averaged eigenvalue density, < g(x) >, where g(x) := 5(x — x u )/N, 

is the microscopic density of the eigenvalues, {x u }^ =1 , of a iV x N Hermitian random matrix. The 
A function f(x) is Holder continuous if 

< c\xi -x 2 \ s , 

for all X\,X2 in the domain of f(x), for a constant c > and < 5 < 1. 
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validity of the mean-field approximation, for large N, is discussed in ||. An easy calculation shows 
that a functional derivative of a w.r.t. v is the density-density correlation function^ : 

* . _ , 5a(x) 5 < g(x) > . . , . . . , . 

(7) C ee (x,t) := -^y = ^ =< g{x) X g(t) > - < g(x)g(t) >, x,t € J 

and must satisfy the obvious sum rules, J ; C ee (x,t)dt = = JjC gg (x,t)dx. Furthermore, from (0) 
C gg (x,t) = C ee (t,x). In this work we explicitly determine C gg (x, t) as a function of the end-points 
of the interval J := U|jti(oj, bj). It turns out that the density-density correlation function can be 
identified with the Bergman kernel of a Riemann surface which is a two-sheeted covering of the 
complex plane. 



2 Determination of the density-density correlation function 

We fix some notation. We consider the hyperelliptic Riemann surface S g of genus g, defined by the 
equation 

S 9 := zeQP 1 , y 2 = JJ( Z - aj )(x - bj) j . 

The projection (y, z) — > z defines S g as a two-sheeted covering of the complex plane C cut along J. 
On S g we define the canonical cycles {a^, /3fc}^ =1 shown in Figure 1. Let 



9-1 

(8) U g {x) ( xP + ^Kjot 

with Kj, j = l,...,g, determined by 

(9 ) r +1[ lM dx = o, j = 1, ...,g. 

Jbj y(x) 

For completeness we first determine A as a function of v. Multiply (||) by \J g (x)/y(x), integrate 
w.r.t. x over J and noting that, 



U fl (x) 



/ ~ yK ~' dx = 1, 

Jj y{x) 



2 Taking a functional derivative of 

< Q(x) >:- f e -{B[ e ]+Jv(x') e (x>)dx')DQ ' 

w.r.t. v(t) gives (Q) 
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we find 



A[v] = -2 / o(t)dt f ^$ln(b g+1 -x)dx+ [ v{x)U9 } x) dx 

J j J j y{x) J j y{x) 



y(x) 



where 



TT(]g(t) 1 



(11) V[J]:= / dt-lnb ff+1 . 

The above integral is along any path connecting and oo + laying on the upper half plane of 
the upper sheet. So in the absence of the external field, A[0]/2 is entirely determined by the end 
points of the conductor. 

Note that 

(12) = 5^y a{x)dx^j = y^(5bja(bj) — Saja(aj)) + J 8a(x)dx = J 5a(x)dx. 

Using similar calculations, 

r a j+i 

(13) / Sa{x)dx = 0, j = 1, . . . , g. 



So performing a variation on (|2|) gives 

9+1 



9+1 , 

(14) <W(x) - 2^ (ln|x - bj\a(bj) - In jx - aj\a(aj)) - 2 J In \x - t\8a(t)dt = 6A. 



Also in this case, we multiply the above relation by D g (x)/y(x), integrate w.r.t. x over J and, by 
(|To|), we obtain 



(15) 5A = [ ^&8v(x)dx. 

Jj y{x) 

Taking a derivative w.r.t. x on (|14|) produces the singular integral equation 

(16) 2P y 



5 a it) , d5v(x) 

v -eft = — x € J. 



x — i cfx 
Now a possible form for 5a (x), x G J, reads, 

(17) Sa(x) = -^--P [ y^^ldt - f [ ^ly(t) ( [ f d ° ) dt. 

y ! K> 2vr2y(x) Jjt-x dt Z^ 2 ^y(x)Jj dt UK ' \J ak (t - s)y(s) J 
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a, 




Figure 1: The cycles {otj, f3j,8j}j =1 . The part of the cycles that lie in the lower sheet are indicated 
by broken lines. 



In order that 8a actually satisfies fllq), cp^ is taken to be a polynomial of degree g — 1; 
g 

(18) ip k (x) = J^jkix 9 ' 1 , k = l,...,g 



with yet undetermined 7^. See ||. With this choice of ipj., ( |T2D and ( |16D are satisfied. Now 
<Pk(x)dx/y(x) is an Abelian differential of the first kind on the Riemann surface S g . If we choose 
jkl in such way that 

(19) / ^ fc X } dx = 8kj, 
J aj y{x) 

then (|l~3l ) are also satisfied. This completes the solution of (|l^). To see that (|l3|) is satisfied, note 
that, 

8a(x)dx = — — / 5a(x)dx, 

bj 2 -/^ 

where 

(20) 8 g = a g . 

Here the cycles 8j, j = 1, ..,g are shown in Figure 1. Integrating by parts in (|i~7| ) and noting that 
y ( aj ) = = 2/(6j) gives, 



(21) 8a(x) =P J C ee (x,t)8v(t)dt, x € J, 
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where 



(22) C M 0M):=g7 o a ^ --V ^^f / ? * V x,te J. 

V 7 eek ^ Si ^27T 2 y(x)(x-i) ^ 2vr 2 y(x) A fc (« - t)y(s) J 

An inspection of (^) shows that sum rules are satisfied. Furthermore, C eg (x,t) can be identified 
with the Bergman kernel of the Riemann surface S g and is therefore symmetric under the exchange 
of x and t, namely C eg (x,t) = C ee (t,x) [g, p218]. 

A more direct way to check the symmetry property of C gg (x, t) is shown below. Let nj(t)dt/y(t) 
be an Abelian differential of the second kind with vanishing a periods 

(23) / ^Mdt = 0, j,k = l,...,g, 

Ja k y{t) 

and with behavior at infinity 

(24) ^Mdt ~ ± ft*'" 1 + 0(t- 2 )) dt, i~oo±. 

y(t) 

The quantity ttj is a polynomial in t of degree g + j: 

(25) 7Tj (t) = r t 9+j + Tit***- 1 + ... + Vjt 9 + a!^- 1 + a 2 t^ 2 + ... + a s , 

where the constants a^'s are determined by ( |23| ) and the constants T^'s determined by ( |24| ) are the 
coefficients of the expansion 



(26) y(z) ~ (Y + V -± + V -± + J) 



z ~ oo 



From the Riemann bi-linear relations the second term in ( |22j ) can be expressed in terms of 7Tj- (t) , 
without involving the constants 7jfe So, 

(27) 2vr 2 C w (x,t) = ^ + ^ + -J-^ ^(2j)r fe _ i 7r i (t). 

y{x)(x-t) y{x){x-ty y(x)y(t) ^ ^ 

The reader can now check that C ee (a;,i) — C gg (t, x) vanishes identically in x and t. For g = 1 this 
symmetry can be established in a straightforward calculation. The kernel, C gg will be used, in a 
later publication, for computing the distribution functions of linear statistics which is of interest in 
random matrix theory. For a discussion concerning linear statistics see S. 
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